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Abstract 

The paper deals with the order statistics and empirical mathematical expectation 
(which is also called the estimate of mathematical expectation in the literature) in 
the case of infinitely increasing random variables. The Kolmogorov concept which 
he used in the theory of complexity and the relationship with thermodynamics which 
was pointed out already by Poincare are considered. 

The mathematical expectation (generalizing the notion of arithmetical mean, 
which is generally equal to infinity for any increasing sequence of random variables) 
is compared with the notion of temperature in thermodynamics by using an analog 
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of nonstandard analysis. 

The relationship with the Van-der-Waals law of corresponding states is shown. 
Some applications of this concept in economics, in internet information network, 
and self-teaching systems are considered. 

Key words: theory of complexity, thermodynamics, Bose gas, Wiener quantiza- 
tion, Maxwell rule, nonstandard analysis, law of corresponding states, distribution 
according to income, Monte Carlo method, queuing theory, World Wide Web. 



2 : 1 Introduction 



Let us recall the history of unrestricted probability theory and quote, in greater length 
than in previous publications, some writings of the great mathematician H. Poincare, who 
was, at the time of these writings, the chairman of the department of probability theory 
and mathematical statistics at the Sorbonne. In his book Science et Hypothese in Chapter 
XI, entitled Calcul des Probability's, Poincare writes: 

"Was probability denned? And can it be defined? And if not, how can we dare argue 
about it? The definition, one may say, is very simple: the probability of an event is the 
ratio of the number of favorable outcomes to the total number of outcomes. 

A simple example shows how incomplete this definition is. Suppose I throw two dice; 
what is the probability that at least one of them will come up with 6? Each dice can come 
up with six different values: the total number of outcomes is 6 x 6 = 36; the number of 
favorable ones is 11, thus the probability is 11/36. 

Such is the correct solution. But I will be equally correct if I say that the throw of 
the two dice can result in (6 x 7)/2 = 21 different outcomes; among these 6 are favorable; 
thus the probability is 6/210 



1 Note that in this (second) calculation the outcome (1,6) is identified with the outcome (6,1), i.e., it 
is counted only once; in other words, we do not distinguish the two dice. We shall see below that, in 
thermodynamics, this corresponds from Boltzman statistics to the statistics of Bose-Einstcin. 
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Why is the first method of calculating the probability better than the second one? 
Whatever the case, our definition doesn't tell us. 

Thus we must complete the definition by saying ". . . to the total number of outcomes 
provided these outcomes are equiprobable" . And so we come to a definition of probability 
in terms of probability itself. 

How can we find out whether or not two events are equiprobable? Isn't that the result 
of some arbitrary convention? If, at the beginning of each problem, we specify what the 
arbitrary convention is, then all will work out well, we will only need to apply the rules 
of arithmetic and algebra and perform the corresponding computation, whose result will 
then leave no room for any doubt. But if we wish to apply this result, we will have to 
prove that our arbitrary convention was valid, and so we return to the difficulty that we 
had tried to circumvent. 

One might say: a bit of common sense will suffice to decide what specific convention 
should be chosen. However, Bertrand, out of pure curiosity, studied the following simple 
problem: 'What is the probability that a chord of the circle is longer than the side of an 
inscribed equilateral triangle?" The great geometer made two conventions, both suggested 
by common sense, and obtained 1/2 in one case and 1/3 in the other. 

The conclusion that apparently follows from the above is that probability theory is a 
useless science, that we should not trust the unclear instinct that we call common sense 
and to which we appeal in order to specify our conventions" (English translation of PQ, 
pp. 115-116). 

The difficulties related to the dependence of probability theory on common sense, 
which were already pointed out by Laplace, were overcome by Kolmogorov in his com- 
plexity theory. We toss a coin but, instead of of talking about the probability of getting 
heads or tails, we perform a sufficiently large number of long trials and then notice that 
most of the possible combinations lie in the interval 1/2 ± y/N, where N is the number 
of trials in the sequences. Here Kolmogorov finds the shortest code that determines an 
individual sequence of trials, and the length of that code will be close to the value of the 
entropy. 

This means, although Kolmogorov does not say it outright, that his approach is based 
on entropy. 

This remarkable achievement, in the case when the values of the random variable are 
not unbounded, completely answers the question asked by Poincare in his philosophical 
work quoted above. 

Here we carry over Kolmogorov' s approach to sequences of random variables in the 
case when the sequences infinitely increase in length. 

First we consider the case of infinitely growing values of the random variable corre- 
sponding to income, i.e., the case of the so-called "Petersburg paradox". A player plays 
on a roulette; he can put his stake on black or red, and begins by putting I dollars on 
black. If the roulette stops on red, he plays black again, increasing his stake to el. Having 
lost, he continues playing black, again multiplying his stake by e, and continues in this 
way until, after m steps, the roulette stops on black. The stake at that point was e m l. 
The player lost 

^ e-1 

k=0 

and won 

e m l. 



His final winnings were therefore 

0.418e m -l, m>2. 

In its time this paradox could not be used by roulette players to obtain garanteed winnings, 
because they did not have sufficient sums of money. But today, when money is electronic 
and high level functionaries in Russia have unlimited credit in state banks, the possibility 
of enrichment in such ways is boundless H. And if society does not have sufficient control 
over the actions of its functionaries, this is very dangerous for the economy. 

If winnings can be regarded as a random variable, more precisely, as the values of some 
random variable £, then the exponential growth of such a random variable is catastrophic. 

In our times, we are nearer to a dangerous economic catastrophe than to the ecological 
one everybody is talking about. In 1567, the (saint) Pope Pius V wrote the bull "De salutis 
gregio dominicis" (On the salvation of the Lord's flock). But neither the inquisition, nor 
the laws of Shiriat, are capable of saving all of society, which consists of people of different 
mentalities and different religious beliefs (see Remark [TU1 below) . 

We would like to point out that historians do not pay sufficient attention to such a 
problem as corruption. Marxist historians regarded ethical wrongdoings and the corrup- 
tion of high level functionaries as too rough an explanation of social problems. 

But already Salustius (86-35 BC) believed that the principal reason for the fall of the 
first oligarchy was the all engulfing corruption. As an example of the implementation 
of strict Roman ethical rules he indicates Caton (on the role of ethical rules in society, 
see [2]). 

Remark 1. In dealing with the order statistics of the income of one of the member 
countries of the European Union, one should consult the statistical data of the country 
before it became a member of the EU. Then one can predict, by using unrestricted prob- 
ability theory, when this country will "fall out into a condensate" like dew condensing 
from the mixture of gases constituting air. 

2 Dependence on the choice of the origin. Notion of 
temperature as mean energy 

The great H.Poincare pointed out the relationship between thermodynamics and temper- 
ature. He wrote: "Let me also mention the kinetic theory of gases and the well known 
hypothesis, according to which it is assumed that each molecule of gas describes a very 
complicated trajectory, whereas the phenomena considered in the mean (the only way 
that they can be observed) are subject, in view of the law of large numbers, to simple 
rules, like the laws of Mariotte and Gay-Lussac. 

All these theories are founded on the law of large numbers, and so the failure of 
probability theory would lead to that of these laws." (English translation of [T], p. 117)o. 

In our work, we construct a probability theory which is intimately related with the 
kinetic theory of gases and is not based on the notion of equiprobability of events. It 
follows Kolmogorov's main idea, which he formulated in the framework of his complexity 
theory. 

2 In France, Germany and other countries there are always ongoing cases in courts of law where such 
wrongdoings by high level functionaries are judged. 

3 Poincare constructed a mathematical model of ideal gas. V.V.Kozlov in [3] proved that the Gibbs 
paradox takes place for this model. Kozlov's proof stimulated investigations of this important problem 
by ] one of the authors (see [I], [5], [5], etc.). 
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Let us clarify the relationship with Bose gas. 

The empirical expectation of a random variable taking the values {Aj} has the form 



oo 

or equivalently 



n * — ' n 

= 1 8=1 



8=1 



M, M = nS, (2) 
n 4 - = n. (3) 



oo 

E 

8=1 

These equations coincide precisely with the equations used to determine the inner energy 
of an ideal gas if A, is the collection of eigenvalues of a certain unbounded self-adjoint oper- 
ator H, known as the Hamiltonian, and N{ is the number of particles at the corresponding 
level: 

oo oo 

J2Ni\i = E, Y, N * = N - ( 4 ) 

8=1 8=1 

The mean energy of one particle is 



-^average 



E 
N' 



As is well known, it is related with the notion of temperature T. 

Indeed, we all know that in a gas the number of particles is huge: in a vessel of volume 
V = 1cm 3 there are N ~ 10 19 particles (at room temperature and under pressure of 
latm). The mean energy of the gas per unit volume is kT, where k is the Boltzman 
constant and T is the temperature of gas in its equilibrium state. 

Therefore, in thermodynamics, a mean quantity corresponding to a huge number of 
tests has been defined in the case when the random variable A assumes unboundedly 
growing values. Its growth will be discussed below. 

Van der Vaals is famous for his empirical relation (the VdV equation). But all sur- 
veys articles of thermodynamics claim that van der Vaals' main discovery is the law of 
corresponding states. This law consists in the following. Real pure9 gases have the same 
temperature T c and the same critical pressure P c . 

D.Mendeleev referred to the critical temperature T c as the temperature of absolute 
boiling. He had in mind that the liquid state is impossible at higher temperatures. What 
is temperature (in the Kelvin scale)? It is counted off from absolute zero. However, the 
absolute zero is unaccessible in the sense that no matter how small the temperature (in the 
Kelvin scale), it still is infinitely far from absolute zero. In this sense it is more natural to 
consider the logarithm of the temperature. Then logT|<r =0 = — oo, and it is obviously seen 
that the distance between the temperature (in the Kelvin scale), no matter how small, 
and absolute zero is an infinite quantity. Thus we are dealing, from the mathematical 
standpoint, with "nonstandard analysis" and an infinite temperature, although the latter 
may equal, say, one hundredth of a degree (in the Kelvin scale). 

If we pass to the relative quantities 

T P 

T r ed 7p~ i Pred p i (5) 



A pure gas is a gas consisting of identical molecules. 
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then it turns out that in these reduced coordinates, the isotherms, isochores, and isobars 
on all the coordinate planes (pressure - volume), (temperature - density) and so on 
behave very similarly for different gases even when the gases have very different critical 
temperatures T c . This is true for nitrogen, argon, carbon dioxide, methane and other 
gases. 

Example 1 . Consider the standard Hilbert space H of functions defined in the closed 
interval [0, tt], and its basis of orthogonal functions 

{e-}n=-oc, (6) 
and the following two self-adjoint operators 

A, = -ih—, A 2 = -2ih—, (7) 
ox ox 

depending on the parameter h — l/k, where k is an integer. 
The initial reference point in the basis {e mx } is n = 1. 
Consider rijp = k and tiq = §. We see that 

A ie in o 1>x « 1, 



^ . (2) 

A 2 e m ° x « 1. 



Obviously, 



e 



i(n—k)x 



Y 

) r 



i n=oo 
-00° 



n=— oo 

and 

I J 71= — OC 

are isomorphic to 

Thus, as h — > 0, the reference point of the infinite basis shifts from the point n — 1 
to the point n = 1/h in the first case and to the point n = l/(2h) in the second 
case A 2 . On the interval [0, tt], the operators A\ and A 2 will act exactly in the same 
way on functions x expressed in the basis ([6]) as the operator A of the usual numbering 
of functions. When we shift the numbering by n — k and n — k/2, we come to rapidly 
oscillating functions. These two spaces are orthogonal to each other and to the initial 
Hilbert space H (in the sense that J sin j^<p(x) dx — > as h — > 0, ip(x) G H). 

Denoting e lkx = ty c and e %k l 2x = ip c , we see that the space of functions corresponding 
to the basis of the form ^/ n /^/ c is isomorphic to the basis of functions ^ n /ip c as h — > 0, 
where ty n = e mx (h = 0, k = oo in the "nonstandard analysis"). 

A similar situation occurs for Dirichlet series and for the eigenf unctions of the 
Schrodinger operator (the Hamilton operator) H) for the oscillator and so oe@ (see [Hj). 

This isomorphism of infinitely large quantities from "nonstandard analysis" explains, 
in a certain sense, the empirical law of corresponding states due to Van der Vaals, who 
discovered a similar "almost isomorphism" . 

It is difficult for a mathematician to explain the ideas of the new thermodynamics 
to a physicist. The most difficult thing to explain is that such a small temperature as 



5 Note that the relationship between the solution of the Cauchy problem for the Schrodinger equation 
of the form e ,Bi $o, where € H is the initial condition, and the solution of the heat equation e~ Ht ^o, 
is described in Section 2 of the paper [7] which deals with Wiener quantization. 
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T = 1/100 K is actually an infinitely large quantity For the expectation ([I]) of a random 
variable taking unboundedly increasing values (as, for instance, in the case Aj = i 2 , 
i = 1,2,...), this fact is quite natural in the case when we decompose the natural number 
M into N summands, as well as in the case of frequency dictionaries [9]. 

Example 2. Let us recall a famous theorem from number theory, namely the solution 
of the "partitio numerorum" problem. In this problem we deal with a natural number M 
that we are to decompose into the sum of iV summands (positive integers), for example 
M = 5, N = 2: 

5 = 1 + 4 = 2 + 3, 

which gives 2 variants of the decomposition: Ai = 2. 

If M = 10 23 , N = 1, there is only one variant: M = 1. If M = 10 23 and N = 10 23 , 
there is also only one variant of the decomposition (the sum of ones): M. — 1. 

Obviously, for a fixed M there exists a number N c for which the number of variants A4 
is maximal (this number is not unique in general). The quantity log 2 M. is called Hartley 
entropy. The point at which it achieves its maximum is where the maximal entropy 
occurs. 

Suppose we are given the decomposition 

M = ai + h a N 

of the number M into N summands. Denote by Nj the number of summands in the 
right-hand side of that equation which are exactly equal to j. 

Then there will be J? . Nj summands, and this sum equals N, since there are N 
summands in all. Further, the sum of all summands that equal j is jNj, since there are 
jNj summands Nj, so that the sum of all the summands is obtained by summing all these 
expressions over j, i.e., it is ^2jjNj, which equals M. Namely 

oo oo 

j2n = n, j2 m * = M - ( 8 ) 

i=l i=l 

The nonuniqueness of the above mentioned maximum and the indeterminacy of their 
quantity permitted Erdos to obtain his result only with precision o(a/M). 

Thus we have the Erdos theorem for the system of two Diophantine equations 

oo oo 

J2^ = N, J^iV 4 = M, (9) 

1=1 8=1 

which asserts that the maximum number of solutions N c of this system is 

N c = $- x M\l 2 log M c + aM l J 2 + o(M c 1 / 2 ), = ir^/2/3, (10) 

where the coefficient a is determined by the formula (3/2 = e~ a/3 ^ 2 . 

If we increase the number N in problem while M remains constant, then the 
number of solutions will decrease. If we count the sums ([§]) from one instead of from zero, 
i.e., 

oo oo 

Y^m = (M-N), ^iV, = iV, (11) 

8=0 8=0 

then the number of solutions will not decrease, and will remain constant. 



Let us explain this fact. The Erdos-Lehner problem [10] consists in decomposing the 
number M c into N < N c summands. 

The decomposition of the number 5 into 2 summands gives two variants. If we include 
the number 0, we obtain three variants: 5 + = 3 + 2 = 4 + 1. Thus the inclusion of zero 
gives us the possibility of saying we are decomposing M into k < n summands. Indeed, 
the decomposition of 5 into three summands includes all the previous variants: 5+0+0, 
3+2+0, 4+1+0 and adds new variants without zero. 

The maximum does not change by much in this case [10] . whereas the number of vari- 
ants does decrease: the zeros will allow the maximum to remain constant, and the entropy 
will not decrease, it will become constant after the maximum is reached. This remarkable 
property of the entropy is what allows us to construct the unrestricted probability theory 
in the general case [IT]. In physics, this effect is identical to the so-called Bose condensate 
effect. 

One may ask: How does the "partitio numerorum" in arithmetic differ from the 
Boltzman- Shannon statistics? If we regard 4 + 1 and 1 + 4 as different variants, then we 
obtain the Boltzman-Shannon statistics. The number of variants AA of the decomposi- 
tion rapidly increases. Thus the "noncommutativity" of addition gives a huge number of 
extra variants of the decomposition, and the Hartley entropy, equal to the logarithm of 
the number of variants, will coincide with the Boltzman-Shannon entropy. 

The entropy arising in arithmetic, because of commutativity and symmetry, is the 
Bose-Einstein entropy. Let us recall the formula for the entropy of a nonequilibrium Bose 
gas, derived in the textbook [12] (formula 54.6): 

S = G A(1 + *i) Ml + nj) ~ % ton,], (12) 

j 

where nj are the average numbers of fillings of the quantum states rij = Nj/Gj. 
We have 

3 3 

£ > x i H -/•>; ( 13 ) 

3 3 

From this, applying the method of Lagrange multipliers to the maximum of the en- 
tropy (fl2|) under conditions ( TT3|) . we obtain^ 



e a+h£ j - 1 v ; 



Here a = —fi/T, b = 1 /T, T is the temperature, \x < is the chemical potential. Thus the 
Lagrange multipliers are expressed in terms of the temperature and the chemical potential 
of the gas. 

Since, according to the Bose-Einstein formulas sums of the form 



OG 

n, a. 

i=0 i=0 



are replaced by integrals, it is impossible to distinguish the point i = itself, the conden- 
sation effect must occur in a neighborhood of this point, not at the point % = itself. But 



6 We repeat the derivation word for word [T2]. A rigorous proof of fractional dimensions was first 
obtained by A.M.Vershik [T3]. Also see [H]. 
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for D < 2, the Euler-Maclaurin formula cannot be applied (and so we cannot pass from 
sums to integrals), and therefore the point i = can indeed be distinguished. 

However, will the extra "particles" condense precisely near the point i = for iV > N c l 
Actually, this question consists in the following: will the number of zeros in the problem 
under consideration as M — > oo be much larger than the number of ones, i.e., will we 
have 

Nt = o(N ) 

even for N > N c 7 

A computer calculation gives a negative answer. Experts in number theory also answer 
negatively. All the more when D > 2 and the Euler-Maclaurin formula can be applied, 
for i 3> 1 we can pass to the intergral and it is impossible to speak of one point % = 0, as 
it is done in the textbook [T2] . 

Actually, together with the convergence to the integral (for instance, in the 3- 
dimensional case) as iV — > oo and — > 0, the bell-shaped function \/Ne~ Nx also 
increases, where 

y^>*>0 (15) 

(cf. nsi). 

In this general argument it is important to understand the meaning of Gj in the general 
case. 

In the mathematical literature, there are two definitions of Gj\ the statistical one, cor- 
responding to quantum statistics, and the one from mathematical physics, corresponding 
to the .D- dimensional Laplace equation in a D-dimensional volume. In exactly the same 
way, the Landau-Lifshits textbook has two definitions of the Bose gas. One appears in the 
part called "Statistical Physics" , the other one (which uses the symmetric eigenfunctions 
of the Schrodinger operator) appears in the part entitled "Quantum Mechanics" . In the 
mathematical literature, the first is named after Weyl, the second, after Courant. Let us 
look at both definitions in more detail. 

The 2.D-dimensional phase space is discretized into a lattice and Gj is determined 
from the formula 

G. = (16) 



Here h is the Planck constant. 

By assigning to Gi the multiplicities of the spectrum of the Schrodinger operator, we 
obtain a correspondence between the symmetric (with respect to permutations) eigenfunc- 
tions of the 
iV-particle Schrodinger equation with the combinatorial calculations of the Bose statistics 
appearing in the book [T2] . 

The single-particle ?/>-function satisfies the free Schrodinger equation with Dirichlet 
conditions on the walls of the vessel. According to the classical Courant formula, the 
number of eigenvalues less than a given A ^> 1 (the spectral density p(A)) has the following 
asymptotics: 

y m D/2 X D/2 

p(A) = rWHlpW (1 + ° (1)) as A^oo (17) 

(V is the D-dimensional volume). 

This is precisely by means of this correspondence that we establish the relation- 
ship between the Bose-Einstein combinatorics [12] with the definition of the A-particle 
Schrodinger equation and the multiplicity of the single-particle Schrodinger equation. 
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The spectrum of the single-particle Schrodinger equation without the interaction po- 
tential coincides up to a factor with the spectrum of the Laplace operator. 

In the case of the Laplace operator and the Schrodinger equation, it makes no sense 
to consider fractional dimensions. In the statistical case, fractional dimensions do make 
sense, but their meaning is different. They determine the average number of degrees of 
freedom, which can be different for molecules of different velocities, so that the average 
value of the degrees of freedom for the whole gas can be fractional. It corresponds to the 
mean energy, which is the temperature. For the unrestricted probability theory that we 
consider here, it is precisely the fractional "dimension" (understood as the average of the 
number of degrees of freedom over the entire "general" population) that makes sense. 

One of the authors proved that both the degeneration energy for the Bose gas and 
the number of particles in the degeneration state coincide up to normalization with the 
critical values N c and E c (i.e., for fi = 0) of the classical noninteracting gas. In the 
-D- dimensional case, the degeneration energy has the form 



E d 

where 



2m 



o e ^~ //Td — 1 



d£ Jtf d Pl ■ ■■dpndVn 

2m (2-kK) d ' v ' 

which implies the coefficients C and A D in the formula 

E d = CA D T 2+7 C(1 + £>/2)r(l + D/2). (20) 
It is more convenient to pass to another normalization 

de' = — — -de (21) 

r(D/2 + i) v ; 



for the energy and, respectively, 



ds " = nm' k (22) 



for the number of particles N. 

Passing to the Van-der-Waals coordinate (jSJ), we obtain the coefficient CA D of the 
form A 7_7c , where 7 C is determined from the experimental values of the quantity Z\^ =0 = 
((jc + 2)/£(7c + 1) and ((x) is the Riemann zeta function. In what follows, we omit the 
tilde on the coefficient A and the index r (reduced) on the parameters T and P. 

After this, it is more convenient to use the polylogarithm Li 2+7 (e _M//T ) for the energy 
and the polylogarithm Li 1+7 (e~ M / T ) for the number of particles: 

£ = T 2+7 Li 2+7 (a), (23) 
iV = T 1+7 Li 1+7 (a), (24) 

where a is the activity (a = e~^^ T ). 

The the compressibility factor is equal to 

v _ E _ Li 2+ » _ M 



NT Lix+ 7 (a) Tn 



The coefficient A is here eliminated. 



3 Basic series and phase transition 



To each degree of freedom we assign its basis sequence, whose terms correspond in ther- 
modynamics to pure gases with identical molecules. Namely, we consider the following 
basis family of random variables {i D / 2 }, corresponding to the number D of degrees of 
freedom. 

The variational series 

6 = * D/2 , (26) 

will be called the basis serie$\ of dimension D. 

Let rii be the value of the outcome corresponding to 

A sample value (estimate) of the mathematical expectation is 

E>«& = g 
n 

We assume that £ — > oo as n — > oo. 

It is useful to introduce the parameter 7 as follows 

D = 2(7 + 1), 7 = y-l- 

We shall consider the one-parameter family of entropies S 7 as the logarithm of the number 
of solutions of the system 

00 00 

^ rtj = n, ^ n & ^ n£ - ( 27 ) 

i=l i=l 

The asymptotic relation n(S) between n and 8 as n — > 00 for which the maximal num- 
ber of solutions of (1271) (up to o(N)) will be called the point of maximum des accumulation 
of the estimate (sample value) of the expectation, while the families of solutions {n^/n 
of equation (j2"7|) will be called typical] the value n for which the maximal deacccumulation 
is achieved is denoted by iV c . 

The binary logarithm of the number of solutions of equation ([27]) is the Hartley en- 
tropy S. 

Thus a, b are the Lagrange multipliers (see (jl2]) - (114p ) for finding the maximum of 
under conditions ( 1271) : 

alS = aaln + b d(£n) , En = M, (28) 

d(£n) = \dS-^dn. (29) 
b b 

Denote the quantity 1/6 by T, and the quantity a/b, by \i (p < 0). In thermodynamics, 
there is an additional term PdV, where P is the pressure and V is the volume. 

Remark 2. In the economic calculation of income (see below), the variable T c will 
be called the minimal wage, while p, = —fi will be the nominal percentage or the desac- 
cumulation coefficient. We also introduce the notation b = 1/T. To the value n = N c 
corresponds T = 1. For a given E, the maximum of desaccumulation is attained at the 
value of n(E), which we denote by iV c . To the volume V corresponds to volume of goods 
Q, and to the pressure P, their price (Irving Fisher's correspondence law). 



7 Similarly to thermodynamics, the pure gases correspond to the periodic table, and each pure gas is 
associated with a certain value Z Cy for example, up to 0.01%. 
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The Hartley entropy S 1 for 7 > has the form 

(30) 

For fi = and for n > N c , the entropy S can be extended continuously by a constant (a 
phenomenon that, in the particular case D = 3, is equivalent to the physical notion of 
"Bose-Einstein condensate" ) . 

The relation n > n c is equivalent to T < 1. Here for n ^> k ^> 1 (see (Tl5!) ). the 
number of values greater than the critical one is X^fc n i = n c- The remaining n — n c 
values accumulate within the interval i < k n which becomes relatively smaller as 
n — )■ 00 and \i — )■ 0. Here k 3> 1 is the support of the bell-shaped function. 

The passage to the integral is accompanied by the appearance of a bell-shaped function 
that progressively approaches the 5-function (Bose-condensate). Therefore in formula fl30|) 
for n > n c or T < 1 the values of dn (and dV in thermodynamics) can be set equal to 
zero. As the result, the entropy of the "liquid" phase acquires the form: 

d(Sn)=TdS. (31) 

On the plane {M, Z = ^}, the "gas" critical isotherm 

M = Li 7c+2 (e-^), Z = P^P, (32) 

Li 7c +i(e ^ ic ) 

is practically the same for all pure gases each of which is associated with Z c = C(7c + 
2)/C(Tc + !)■ F° r argon Z c = 0.29, for methane Z c = 0.29, for oxygen and other gases, 
the value of Z c is determined experimentally. The "base" spinodal (i.e., the set of points 
for fi = 0) on the plane {M, Z} is the interval Z = Z c , M = {0, 1}. The isotherm T < 1 
demonstrates a good coincidence on the interval —00 < /x < //; up to the phase transition 
"gas-liquid" for \i = 

On the "liquid" isotherm, the value n = T 7c+1 £(7 c + 1) is constant. Therefore, it is a 
straight line passing through the point M = T 7c+2 , Z = Z c (the origin is M = 0, Z = 0). 

Remark 3. The Bose condensation in income statistics means that when we increase 
the number of trials n > no, we necessarily come to the low income group if the general 
population is the country's population, while the order statistics is determined by the 
value of the random variable (expressing the income of each citizen of the country; see 
Sec. ED . 

4 Negative values of the parameter 7 

Under the assumption that the transition to the liquid phase does not occur at T = 1, we 
equate the chemical potentials \x and jx for the "liquid" and "gas" phases on the isotherm 
T = 1. 

Then we determine the point fi, i.e. the point of transition to the "liquid" phase for 
T < 1, by equating the chemical potentials of the "liquid" and "gas" phases. 

Remark 4. Consider the case 7 < and jx = 0. 
Let us introduce the notation 




S 



n 



Z( 7 )(2 + 7 ) + 



T 



1 1 



For 7 < and fi = we have (according to ([I 



OO ( oo oo ^ 



& = (33) 



Since the function f(x) = x 1 F{bx) monotonically decreases, we have 

00 00 /*00 /*OG POO 

J2fF(bj) = J2f(j)< f(x)dx = x 7 F(6x) rfx = r 7 " 1 / x^F(x)dx (34) 

(the Nazaikinsky inequality). 
Thus, we have 

iV| A=0 = A^6- 1 C(l-7) + O(6- 1 ^), 7 <0, 6 = ^, (35) 

where £ is the Riemann zeta function. 

In this section we find the liquid isotherm-isobar point as the quantity x = —fi/T 
slowly tends to zero. 

First, we recall that although N c is large, it is still finite and hence we must use the 
parastatistical correction. We recall the relation 

^ /1-exp^A^ 



for the fl-potential. 

We use the Euler-Maclaurin formula with the parameter 7 taken into account and 
obtain 

y(—>: k Ji\ = i _r (_}: *_w+i? 

\g6j+x gbkj+x J Oi Jq \gbx+x Y gbkx+y -y / ' 

where a — 7 + 1, k = N and 6 = 1/T. Here the remainder i? satisfies the estimate 

1 r°° \ h 
\R\<~ \f(x)\dx a , where f(x) = — ^r— . 

We calculate the derivative and obtain 

bk 2 e k( - bx+>() be hx+K 
■> ^ X > ~ ^ e k(bx+x) _ rj2 _ ^bx+y _ rj2 ' 



1*1 <i 



dy a . 



We also have 



+ ip(y), where ip(y) is smooth and \ip{y)\ < C(l+ 2 - 



( e i/ _ 1)2 y 2 
Substituting this formula into (|37|) . we obtain 



R\<- \k 2 ifj(k{y + x))-tfj{y + x)\dy c 
a Jo 



< k 



2-a 



ky 



Ck 



2-a 



\i>(y)\dy a + J ^(y^dyKj— 



1? 



with a certain constant C. For example, it x ~ (In A;) -1 / 4 , then \R\ has the estimate 
\R\ ~ 0((ln k)( 2 ~ a ^ 4: ). Therefore, we have k = N c and T = T c and obtain the following 
formula for the integral at \x = 0: 

A7-7c f £d£ a A 7 " 7c y 00 7]dr] c 



ar( 7 + 2) J e fe « - 1 J & - 1 ' 

where a = 7 + 1. Thus, 

1 / A 7 ~ 7c /"°° £c/£ Q ^ 



IVW V" r (7 + 2 ) 7o e« - 1 

We obtain (see [T7] ) 



e b§ _ 1 e feb? _ 1 j * 5a y o y e t _ 1 £ 



1 r (-- 1 \ _ z 00 r _^ k a \ 



■ c(7 V-°-i)- 



b° 

Since 3> 1, we finally have by setting k = N\^/ T=0 ^, 



iV| A /T=o(i) = (A 7 " 7c c( 7 )) 1/(1+7) T, where c( 7 ) = / (7-^— y)^^- (3 



The relation AT = T 7c+1 £( 7c + 1) agrees with the linear relation A" = A(j)T, where 
A( 7 ) = (A 7 ->c( 7 )) 1 /( 1+7 ), for P < 0. 

We can normalize the activity a at the point T c and can find do, by adjusting the 
liquid and gas branches at T c for the pressure (i.e., for M) so that the phase transition 
on the critical isotherm could not occur at T = 1. 

Further, we normalize the activity for T < 1 by the value of calculated below. 
Then the chemical potentials (the Gibbs thermodynamical potential for the liquid and 
gas branches in thermodynamics) coincide, and hence no "gas-liquid" phase transition 
occurs at T = 1 (cf. Sec. 5). 

Now to realize the isochor-isotherm of "incompressible liquid" , we must also construct 
it a density lesser than N c = C( 7 c + 1), namely, for 

N(T) = T> +1 C( 7c +l). 

We obtain the value 7 (T) from the implicit equation 

A( 7 )=r*C(7c + l)- 

We thereby obtain the spinodal curve (i.e., the points where \x ~ T(ln A^(T)) -1 / 4 ) for 
each T < 1 in the region of negative 7: 

A -(7-7c)/(i+7) c ( 7 ) V(l+7) = T 7c C(7c + 1). (39) 

From the two values of 7 of the solution (|39l) we take the least (the largest in the 
absolute value) and denote it by 7 (T). In particular, for T = 1 we denote 70 = 7 (1). 
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The second solution can be realized for some metals (mercury, cesium, etc.) and 
corresponds to the metastable state. The actual interaction between particles gives a 
more significant deviation of the second solution from 70 (see [T6]). 

Let a g = e~^l T be the gas activity, and let a z = e~^l T be the liquid activity. We 
give the condition that M and the activities coincide at the point of phase transition (see 
Sec. 5): 

T* Li 2+7c (a 9 ) = A-M^T-I^)I Li 2 _ |7(T)| (f) , (40) 

\ a o/ 

Li 2+7o (a ) = 1, (41) 



C(2 + 7o) 

a g = ^. (42) 
a 

These two equalities determine the value of the chemical potential /i = pL = T In a g at 
which the "gas-liquid" phase transition occurs. 

Let T = min_ 1<7<0 ^(7). Thus, for each T < T < T c we obtain the value of the 
liquid reduced activity a r = ai/a (where ai is the liquid activity), which corresponds to 
the Van-der-Waals normalization (see (jSJ))- 

In thermodynamics, the critical values T c , P c , and p c were experimentally calculated 
almost for all gases, and hence the critical number of degrees of freedom can be given in 
advance. The parameter A (1.6 < A < 3, T > 1/3) is obtained from the coincidence of 
binodals at the triple point determined experimentally. 

Remark 5. In the probability theory and economics, these quantities can be deter- 
mined only by indirect methods (for example, by examining the income growth rate (see 
below) or by an analysis of a sufficiently large set of debt-duration pairs (duration is the 
period of debt repayment); see [IE])- But the analogy with the Van-der-Waals law of 
corresponding states allows us to hope that the scenarios of "phase transitions" must be 
sufficiently similar for different states (see [1$]). 



5 Infinitely small dissipation taken into account and 
Maxwell rule 

A.S.Mishchenko and one of the authors [20] considered the case of transition to the two- 
dimensional Lagrangian manifold in the four- dimensional phase space, where the pres- 
sure P and the temperature T (intensive variables) play the role of coordinates and the 
extensive variables, the volume V and the entropy S play the role of momenta of the La- 
grangian manifold, where the entropy is the action generating the Lagrangian structure. 

It turns out that this complicated transformation, which leads to the Gibbs thermody- 
namical potential, can be performed only by using the extension to the region of negative 
values of 7. Then one can justify the Maxwell transition by introducing a small dissipa- 
tion (viscosity). The introduction of infinitely small dissipation allows one simultaneously 
to solve the problems of critical exponents without using the scaling hypothesis, which 
underlies the renormalization group method. Let us show this. 

There is no viscosity in thermodynamics. However, in general, thermodynamical equi- 
librium cannot occur without infinitely small dissipation. Therefore, it is natural to 
introduce a small viscosity and then to let it tend to zero. 



14 



The geometric quantization of the Lagrangian manifold (see [21], §11.4) is usually 
associated with the introduction of the constant H. The case of pure imaginary number h 
was called the Wiener (or tunnel) quantization by one of the authors [22] . 

Let us apply the Wiener quantization to thermodynamics. The thermodynamical 
potential G = fiN is the action § = j p dq on the two-dimensional Lagrangian manifold A 2 
in the four-dimensional phase space q x , q 2 , Pi, P2, where qi and q 2 are the pressure P and 
the temperature T, respectively, p\ is equal to the volume V, and p 2 is equal to the 
entropy S with opposite sign. All other potentials, i.e., the internal energy E, the free 
energy F, and the enthalpy W are the results of projection of the Lagrangian manifolds 
on the coordinate planes pi, p 2 



W 



E = - J qdp, q = {qi,q2},P={pi,P2}, 

- J(q 2 dp 2 + q 1 dp 1 ), F = J(q 1 dp 1 -q 2 dp 2 ). (43) 



In the Wiener quantization, we have 

N — e ^ V — e ^ S - e ® 
<9// dp' &T 

Hence the role of time t is played by ln(— n/T) 

H d d 



G = /jlN ~ e 



Td{fi/T) <91n(/i/T)' 



Note that the tunnel quantization of the Van-der-Waals equation (VdW) gives the 
Maxwell rule as e — > (see below). 

As is seen below, the critical point and spinodal points are focal point and hence do not 
"pass" into classics, i.e., into the VdW model, as e — > 0. The spinodal points, similarly to 
the turning points in quantum mechanics, are approximated by the Airy function, and the 
critical point, which is the point of origination of two turning points (two Airy functions) 
is approximated by the Weber function (see [23]). The point of shock wave origination in 
the Burgers equation is expressed precisely in terms of the Weber function as e — > 0. If 
the transition as e — > is performed outside these points, then the VdW-Maxwell model 
is obtained. But the limit transition is violated at these points themselves. Therefore, 
the so-called "classical" critical Landau exponents [12J turn out to disagree seriously with 
experiments. The Weber function gives singularities of the form e~ 1//4 , and the Airy 
function, of the form e~ 1//6 . 

We now consider the Burgers equations in more detail. 

Consider the heat equation 

du e d 2 u m , , 

a=2a?- ieR - i£0 - < 44 > 

where e > is a small parameter. It is known that all linear combinations 

u = Ai^i + X 2 u 2 (45) 

of the solutions Ui and u 2 to Eq. fT44|) are solutions of this equation. 
We perform the change 

u = exp(—w(x,t)/e) (46) 
1F> 



and obtain the nonlinear equation 

dw 1 f dw\ 2 e d 2 w , , 

ST + 2UJ -2l? = - (47) 

which is called the integrated Burgers equation H. Obviously, any solution Ui of Eq. (jHJ) 
is associated with the solution Wi = — elnui of Eq. ( I47p . i = 1,2. The solution (T4"5l) of 
Eq. (jHJ) is associated with the solution of Eq. f|47|) 

w = — e In (e e + e £ j , 
where /Xj = — elnAj, (2 = 1.2). Since 

lim w = minfiyi, W7 2 ), 

we obtain the (min, +) algebra of tropical mathematics |24j . 

To determine the solutions of t > i cr , Hopf proposed to consider the Burgers equation 

dv dv e d 2 v . . . ,, . 

¥ + ^-2^ = °- «U =«.(«)■ (48) 

The (generalized) solution of the equation 

% +P % = {) ' Pl*=o=Po(x) (49) 

is the function p gen = lim £ ^. v (the Riemann waves) . 

The solution v of the Burgers equation is expressed in terms of the logarithmic deriva- 
tive 

d 

v = —s— In u (50) 
ox 

of the solution u of the heat equation 

¥ = l£' = exp { - \ L Mx) dx \ (5i) 

Thus, the initial problem is reduced to studying the logarithmic limit of the solution of 
the heat equation. It is known that the solution of problem (15 ip has the form 



« = (2net)- 1 / 2 J exp|-^(x-0 2 + 2t^ p (£) <%) / 2th} d£. 



(52) 

The asymptotic behavior of the integral (152p is calculated by the Laplace method. For 
t < t CT , we have 



u = (|J|- 1/2 (£(x,t),t) + 0(e)) exp | - ^(x,t)|. (53) 



Here S(x,t) = f_ pdx, and the integral is calculated along the lagrangian curve A*, 
and the point r(x) is a point on A*. For t > t CI , there are three points r\(x), 7*2(2;), 7*3(2;) 



8 The usual Burgers can be derived from Eq. (|47j) by differentiating with respect to x and performing 
the change v = dw/dx. 
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on A* whose projections on the axis x coincide or, in other words, the equation Q(t,£) = x 
for x G (xx,X2) has three solutions £±(x,t), &(x, t), ^(x,t). 

We introduce the notation S(x,t) = f^pdx for x < Xi, x > x 2 , and S(x,t) = 

min(5'i, S2, S3), Sj = fj^ p dx, where j = 1, 2, 3 for x G [x%, x^]. 

The above arguments allow us to obtain a generalized discontinuous solution of prob- 
lem P9|) for times t > t cr . It is given by the function p = p(x,t) determining essential 
domains [22J of the curve A*. We note that, in particular, this implies the rule of equal 
areas known in hydrodynamics to define the front of the shock wave whose evolution is 
described by Eq. ( H9|) . We note that this exactly corresponds to the Maxwell rule for the 
VdW equation. 

The solution v = v(x,e) of the Burgers equation at the critical point x = p 3 is 
calculated by the formula 

mnu( X ) g exp{ f g i g%g 
( '> dx J~exp{=^R' 1 ' 

As x — > 0, after the change -h^ = rj, we obtain 

v(e, x) -$ x ^q \fe ■ const. (55) 

What does this mean from the standpoint of classics and classical measurements, i.e., 
in the case where the condition called a condition for being "semiclassical" in the book [25] 
is satisfied (i.e., in the case where we are outside the focal point)? For the Laplace trans- 
formation this means that we are in the domain where the asymptotic Laplace method 
can be applied, i.e., where 

1 f°° px-S(p) 

u(x) = —}= I e z dp. (56) 
Moreover, the solution of the relation 







X = (57) 

is not degenerate, i.e., 7^ at the point || = x, and the reduced integral (156"]) is 

bounded as e — > 0. For it to have a zero of the order of e 1 / 4 , we must integrate is over x 
with the fractional derivative D~ l / A . The value of D~ l l A applied to unity, D~ l l A l, is 
approximately equal to x 1//4 . 

According to the indeterminancy principle [7], the correspondence between the dif- 
ferentiation operator and a small parameter of the form D — > 1/e is preserved for the 
ratio — £ du /J >x ; although the leading term of the asymptotic expansion in the difference of 

d 2 u/dx 2 , (du/dx) 2 ■ ,. • , J 

— '- — and - — 2 is eliminated. 

u u 

In the case of thermodynamics, the role of x is played by the pressure P and the role 
of the momentum p is played by the volume V. Therefore, V ~ P 1//4 , i.e., 

P c ~ (V - V c )\ (58) 

This is precisely the jump of the critical exponent. The other critical exponents can be 
obtained similarly [7J. The comparison with experimental data is discussed in the same 
paper. 
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6 Independent events, conditional probability, and 
logarithmic scaling 

Definition 1 . We will call two events independent if the Hartley entropies correspond- 
ing to the events add together. 

Let us consider the conditional probability for independent events. 

It is well known that the conditional probability of an event under a condition is the 
probability of the simultaneous occurrence of the event and the condition normed by the 
probability of the condition: 

where A is the event and B is the condition. 
The formula 

k k 

P{B) = P{B P| AA = £ P(£|A)P(A), J] A t = Q (59) 

i=l i=l 

for nonintersecting events {A{\ is called the total probability formula. 

A similar formula holds in the new conception of probability. It suffices to prove it 
for k = 2. This, according to Kolmogorov's complexity theory, must be the equality of 
entropies with conditional probability {a, (3 = 1 — a} taken into account. Namely, for 
jl = we have 

Z c ( lc + 2) = a( 7l c + 2)Z{ + /3( 72 c + 2)Z C 2 . (60) 

where Z c = ((7 + 2)/C(7 + I), a + (3 = 1. 

From these formulas, we obtain the value of 7c for "mixtures" such as those given 
by ( l5"9~j) by using the basis variational series. 

The above construction can easily be generalized to the general case of formula ( 159]) . 

Logarithmic scaling determines, in the order statistics (variational series), a bijection 
between the numbers A > and In A = a, B > — > b = InB. To the product AB 
corresponds, under logarithmic scaling, the sum In A + In B. 

Addition under logarithmic scaling. 

The summing algorithm is given by 

A®B = — ln( n « +nf ) (61) 
Inn 

with the restriction (1601) on the growth of n taken into account. 

By means of this bijection, we can encode any numbers A and B, as well as the 
semifield AB and the correspondence A + B — > A © B. 

The addition flBTj) as n — )■ 00 undergoes a "phase jump" from the "gas" state to the 
"liquid" state, i.e., from an "open state" to a "closed" one. 

Remark 6. Note that the volume V, one of the main notions in thermodynamics, 
is absent from our considerations, and the role of density is played by the number of 
particles N. During the phase transition gas-liquid, the density undergoes a jump and 
therefore, in our case, so does the quantity N. In thermodynamics, it is natural to assume 
that the number of particles is constant, and the jump in density is due to a jump in the 
volume V. Then the chemical potential fi, multiplied by the number of particles, is equal 
to the Gibbs thermodynamical potential, which remains continuous at the point of phase 
transition (the Maxwell rule). Therefore, the chemical potential of the liquid and the gas 
phases must coincide. 
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Remark 7. This entropy, according to Kolmogorov's complexity theory, is intimately 
linked with encodings of minimal length. Thus the phase transition of the "gas" en- 
tropy ( |30j) into the "liquid" entropy ( 13T|) is the transition to a new encoding. We cannot 
decrease the number of trials: the number of trials is irreversible. In our case an increase 
in the number of trials in the state of "liquid phase" gives us nothing, because the extra 
trials fall out into the "Bose condensate" of the "Bose liquid" . 

Remark 8. In thermodynamics, the volume V is also introduced, but the geometry of 
the volume V is neglected in principle, i.e., the problem is invariant under the variations of 
its shape. There is no external field in this problem, and hence the one-particle potential 
in the Gibbs distribution must have the form (p(r 3 /V) if the molecule radius is "infinitely 
small" . And since the volume is large, the leading term as V — > oo has the form 

ar 3 , . 

- — ■ < 62 > 

where a is a parameter which can be expressed in terms of a new parameter q. The 
parameter q is called the heat of transition from the gas phase into the liquid phase [12] 
as T — > 0. Since the additional parameter must be independent of the shape of the vessel, 
it suffices to take the volume of a ball of radius R for determining this parameter. 

Taking account of this potential changes the pressure and the compressibility factor Z, 
which allows one to take account of the reflection from the vessel walls in the general form. 
We have 

v ' S{i/B * 3 - whs r" e ~" d « = whs ^' mk ~ ■ (63) 

where £ = q{fl- & = 1/T, f3 c = 1/T C . 

Hence the kth term in the definition of polylogarithm 

CO u 



becomes 



whence 



fc=i 

1 (e« fc - l)a k 
1 



P = T~ (Li 3+7c (ae-«) - Li 3+7c (a)} , (64) 



(65) 



t 

z _ Li 3 + 7e (ae~ g ) - Li 3+7e (q) 
Li 2+7c (ae-«) - Li 2+7c (a)' 

Now the condiiton for the spinodal in the domain 7 < becomes 

m = rpjc I Li2+ 7c (e~ € )^— Li2+ 7e (l)| ^ (g6) 

7 Applications 

7.1 Distribution according to incomes 

One of the most important problems of contemporary mathematical statistics is that the 
"general population" of measurable random variables is very large, and the samples are 
also large, but much smaller than the general population. 
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In the USSR the sample used in social questions was usually 3.5 thousand people. 

As an example of the general population, one could consider the set of all citizens of 
the Russian Federation. We would like to calculate their mean income. 

There is a quantity fixed by the state: the minimal wage and the nominal percentage jx. 

On the other hand, the income of many citizens is less than the minimum wage, e.g., 
the family members of the worker (dependents), as well as the unemployed. If we construct 
the order series, the "order statistics" according to increasing incomeqj and slit it up into 
clusters of 1-3%, then, denoting the values of the random variable of ^-income by we 
obtain the variational series 

'k 

Suppose that the number of trials (size of the sample) is n, while n« is the number of 
outcomes for which the number &/T is obtained; then, according to the usual definition, 
the estimate of the expectation of the sample, multiplied by n trials, is equal to 



Thus the estimate of the expectation of the sample is 



E = K ±. (68) 



n 



Obviously, if the £j take unboundedly large values, then a randomly chosen sample 
cannot serve even as an approximate "estimate" of the correct expectation. 

Thus, in determining the average per capita income, we must divide the general pop- 
ulation into parts, distinguishing the region that characterizes (in the mean) the studied 
stratum of society, for which we accept the hypothesis of homogeneous distribution. Some- 
times one considers separate households in their development (longitudinal studies). 

Sometimes one also uses the following methodology for measuring the distribution 
of the population according to income. Here the initial data are regular reports, not 
supported by any documents, of 50 000 people, of which 2000 are muscovites. The choice 
of these respondents is purposely left unchanged (new ones are added only when some 
of the previous ones disappear). There aren't any very rich people among them, because 
participation is purely voluntary. There may be some people among them who live below 
the poverty line, because the reports are modestly paid for. 

The obtained data is extrapolated by comparing the socio-economic characteristics of 
the respondents (region and type of city town or village, gender, age, size of the family, 
profession, etc.) are compared with the data given by census. 

No estimates of the precision of such studies are made, either in Russia or abroad. 
There is only the comparison with "theoretically expected" distributions (lognormal with 
a Pareto tail). Pareto was the first to propose the power distribution as the description 
of the distribution of incomes in accordance with empirical data. There have also been 
theoretical discussions about the possibility of interpreting the observed discrepancies by 
the presence of numerous maxima. 

However, it is impossible to partition an infinite set into a finite family of finite sets. 

Therefore, in our case it is natural to call the estimate of the expectation the empirical 
expectation. 



9 This can actually be done only in the countries where the incomes can be calculated from the paid 
taxes, which means that the incomes are not hidden as a rule. 
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We have called the quantity 

lim = D (69) 

n->oo Inn 

the dimension in the situation corresponding to quantum mechanics and the number of 
degrees of freedom in the statistical case. 

It turns out that for any M (under our assumption (1691) there exists a critical number 
of measurements N c M such that, when the number of measurements increases above 
this critical value N > N c with our asymptotic precision (e.g., 0(N a ), where a < 1), we 
will reach the level of the unemployed or very low income groups, which within the limits 
of our accuracy does not influence the value of N c . 

Remark 9. Let us explain what we mean by a closed system. This is not the analog 
of the situation the ruble in the USSR, although the ruble was not convertible. The USSR 
did not have a market economy, in other words, it did not have an economic system based 
on entropy. 

Another type of closed system is the closed system of currency in postwar Great 
Britain. At that time the exchange of the pound for other currencies was forbidden. 
Similarly, in Israel, the exchange of the shekel was only allowed within the limit of$500. 

A rough analogy with liquid and gas may be carried out here. If we open the cork 
of a bottle of gas, the gas will leave the bottle, unlike the contents of a bottle of water. 
The surface tension, which is the boundary of the liquid, is similar to the boundary of a 
country with a closed system of currency. 



7.2 Problems of queuing theory on the Web 

Here for the general population we take all the communicators: computers, i-pads, i- 
phones, smart cell phones, etc. 

When there is only one communication channel, then a message sent before the others 
will arrive first. Here the order is preserved, and the Shannon theory is applicable, in 
particular we have the Shannon entropy, which is equal to the Boltzman entropy. But 
if the messages are gathered in packets, then this principle can no longer be applied, 
and we come to a specific situation. In problems of queuing theory, the value of the 
random variable defined as the number of dispatched packets, corresponding to the given 
general population, grows without bound. We have in mind the order statistics and the 
corresponding variational series (see [TT]). 

A supercomputer connected to the World Wide Web is similar to a member state in 
a world wide union of states. Sufficient conditions for the computer's dropping out of the 
Web follow from the unrestricted probability theory described above. 

Remark 10. Under the total reign of the World Wide Web, when even high level 
functionaries must confess to it under control of its lie detector, ethical rules prescribed 
by various religions as "rules of the game" and ensuring the existence of human society, 
can be monitored via the Web. 



7.3 Relationship with the Monte-Carlo method 

As an example of the relationship between the Monte-Carlo method, let us describe a 
specific problem for a self-teaching audio-guide, a device that describes paintings in art 
museums. This device is oriented to three degrees of freedom: motion in the direction 
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"North-South", motion in the direction "East- West" and rotation about its axis. How- 
ever, the device must also take into account the directions along which there are no 
paintings (e.g., the directions of certain corridors). 

On the one hand, three degrees of freedom give values of the random variable that 
grow like {i 3 ^ 2 }. On the other hand, the monitoring of corridors without paintings is 
naturally ensured by a separate program in its self-teaching software. 

The application of the Monte-Carlo method for choosing a pseudorandom variable in 
the computer language Java preserves the number of degrees of freedom. For a sufficiently 
large number of trials (i.e., visitors to the museum), the device's software, supposedly 
adapted to the picture gallery, can become, as the result of self-teaching, a software 
adapted to the architecture of that specific museum. And this is the passage from an 
open system to a closed one. 

If the usual rounding procedure is used, i.e., if pseudorandom variables of some finite 
length are considered, then the testing of the initial program is not sufficiently perfect, 
and hence the device cannot meet all requirements of the consumer. In the example of 
an art museum, this fact is not catastrophic. Only the reputation of the firm producing 
this device will suffer. But if such a device is used as a guide of a blind man, then such a 
rounding imperfection can lead to a catastrophe. 
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